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ABSTRACT 

In a universe with a cosmological constant, the large-scale gravitational potential 
varies in time and this is, in principle, observable. Using an N-body simulation of 
a ACDM universe, we show that linear theory is not sufficiently accurate to predict 
the power spectrum of the time derivative, $, needed to compute the imprint of 
large-scale structure on the cosmic microwave background (CMB). The linear part 
of the $ power spectrum (the integrated Sachs- Wolfe effect or ISW) drops quickly 
as the relative importance of Q\ diminishes at high redshift, while the non- linear 
part (the Rees-Sciama effect or RS) evolves more slowly with redshift. Therefore, the 
deviation of the total power spectrum from linear theory occurs at larger scales at 
higher redshifts. The deviation occurs at k ~ 0.1 h Mpc -1 at z = 0. The cross- 
correlation power spectrum of the density 5 with $ behaves differently to the power 
spectrum of 4>. Firstly, the deviation from linear theory occurs at smaller scales (k ~ 1 
h Mpc -1 at z = 0). Secondly, the correlation becomes negative when the non- linear 
effect dominates. For the cross-correlation power spectrum of galaxy samples with the 
CMB, the non-linear effect becomes significant at I ~ 500 and rapidly makes the cross 
power spectrum negative. For high redshift samples, the cross-correlation is expected 
to be suppressed by 5 — 10% on arcminute scales. The RS effect makes a negligible 
contribution to the large-scale ISW cross-correlation measurement. However, on arc- 
minute scales it will contaminate the expected cross-correlation signal induced by the 
Sunyaev-Zel'dovich effect. 



> . 1 INTRODUCTION 



The most intriguing topic in contemporary cosmology is the 
nature of the dark energy which appears to dominate the 
energy density of the Universe at late times. Strong evi- 
dence for the existence of dark energy comes from both the 
combined analysis of the cosmic microwave background ra- 
diati on (CMB) and the ga l axy large-scale stru cture (LSS) 
(e.g. lEfstathiou et al H2OO2I ; ISpergel et afll2003l), and from 
high redshift type la supernovae (e.g. Riess et"all 1 19981; 
iPerlmutter et all 1 19991) . Both of these techniques infer the 
presence of dark energy from geometrical measures. A com- 
plementary probe of dark energy is provided by techniques 
that measure the dynamical effect of dark energy through 
its influence on the rate of growth of structure. Large deep 
galaxy redshift surveys (like the EUCLID, the ESA Mis- 
sion to Map the Dark Universe, and the JDEM, the Joint 
Dark Energy Mission) are being planned that will exploit 
the redshift space anisotropy of galaxy clustering, caused by 
coherent flows into overdense regions and outflows from un- 
derdense regions, to measure directly the growth rate as a 
function of redshift. 

The Integrate d Sachs-Wolfe (ISW) effect 
jSachs fc Wolfe! Il967l ). in which the decay of the large- 
scale potential fluctuations induces CMB temperature 



perturbations, provides another measure of the dynamical 
effect of dark energy. In principle, the ISW effect could 
be detected directly in the CMB power spectrum at very 
low multiples. In the ACDM cosmology, it would boost 
the plateau in the power spectrum at I ~ 10. However, 
as the increase of the power is not large in comparison to 
the cosmic variance, it can not be unambiguousl y detected 
even in the WMAP data l|Hinshaw et all 120081 ). A more 
sensitive technique is to search for the ISW signal in the 
cross-correlation of the LSS with the CMB. As the expected 
signal is weak and occurs on large scales, a very large galaxy 
survey is needed to trace the LSS. Currently individual 
detections based on surveys such as APM, 2MASS, NVSS 
and SPSS (e.g. iFosalba et all |2003|; lAfshordi et all 12004 
Fosalba fc Gaztanaaa 2004: Padmanabh an et all 12005 



Cabre et al. 



2006; Mc Ewen et ail 120071 ; iRassat et al 



2007 



Raccanelli et al. 2008) a re not of ve ry high statistical signif- 



icance (but see Granett et al.ll2008h . There are also analyses 
of the ISW cross-correlation, which combine multiple galaxy 
survey samples and achieve ~ 4cr detection of the ISW 
effect (Ho et al. 2008; Giannantonio et al. 2008). These 
measurements may be the best that can be obtained before 
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the next generation of surveys (BOSS^, Pan-STARRS]0) 
come to fruition and make redshift tomography possible. If 
such surveys are to place robust, meaningful constraints on 
the properties of the dark energy it is important to take full 
account of other processes beyond the (linear) ISW effect 
that may contribute to the cross-correlation signal. Here, 
we focus on deviations caused b y non-linear gravitati onal 
evolution, the Rees-Sciama effect l|Rees fc Sciamalll968T ). 

Other processes are known to contribute to the cross- 
correlation signal. First, the thermal Sunyaev-Zel'dovich 
(SZ) effect (|Sunvaev fc Zeldovich| [l972) caused by hot ion- 
ized gas in galaxy clusters induces an anti-cross-correlation 
signal which can cancel the ISW effect on small scales. 
Its statistical contribution can be modelled and subtracted 
given the value of as (the rms linear mass fluctua- 
tions within a sphere of 8 ft -1 Mpc) which determines 
the abundance of galaxy clusters (e.g. | White et ail 1 19931 ; 
iFan fc Chiuehll200ll ; iMei fc Bartle"ttll2004 ). Also, since the 
thermal SZ effect is frequency dependent, it can be sub- 
tracted in frequency space given sufficient spectral cover- 
age. Second, the redshift dependence of galaxy bias, if not 
properly taken into account, can introduce systematic ef- 
fects in the determination of dark energy parameters. Other 
effects such as lensing magnification and the Doppler red- 
shift effect can also boost the cross-c orrelation signal, but 
are only important at high reds hift l|Loverde et al.l 120071; 
iGiannantonio fc Crittenden! [2007). These effects are well 
documented and can be calibrated and removed. 

In this paper we will solely explore the contribution 
of the non-linear terms, or the Rees-Sciama (RS) effect, on 
the cross-correlation signal. The RS effect arises from th e 
non-linear evolution of the potential (|Rees fc Sciamall 19681 ). 
It is belie ved to be much smaller than the CMB signal at 
all scales (jSeliaklll996l ; iPuchades et al.ll2006h . Indeed, com- 
pared with the CMB power spectrum, the RS effect is orders 
of magnitude lower. Also, compared with the complete inte- 
grated ISW power spectrum, the R S effect has bee n shown, 
using the halo model approach (e.g. ICooravll2002al lbl). to be 
unimportant at I < 100. However, the RS effect has not been 
taken into account in cross-correlation analyses and it is im- 
portant to assess its importance ahead of the completion of 
the next generation of large deep galaxy surveys. 

We use a large N-body simulation to investigate the ef- 
fect of the non-linear contribution on the interpretation of 
the ISW cross-correlation signal. We use the 488 3 - p article 
L-BASICC simulation described by lAngulo et all (|2008l ) 
which, with a box size of 1340 ft -1 Mpc, is ideal for this 
purpose because not only does it enable us to extrapolate 
our analysis to non-linear scales at different redshifts, but it 
includes the very large scale power necessary to check the 
agreement with linear theory. The cosmology adopted in the 
L-BASICC simulation is ACDM, with Q A = 0.75, fi m = 
0.25, n b = 0.024, a s = 0.9 and H = 73 km s _1 Mpc -1 . 

The paper is organised as follows. In §2, we compute 
the power spectrum of the ISW plus RS effects from our 
simulation and compare them with linear theory. In §3, we 
analyse these two effects in terms of the cross-correlation of 



1 http:/ /www. sdss3.org/cosmology.php 
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the LSS with the CMB. Finally, in §4, we discuss our results 
and present our conclusions. 



2 TIME DERIVATIVE OF THE POTENTIAL 

The integrated Sachs-Wolfe effect results from the late 
time decay of gravitational potential fluctuations. The net 
blueshift or redshift of the CMB photons caused by the 
change in the potential during the passage of the photons 
induces net temperature fluctuations of the black body spec- 
trum, 
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where $ is the time derivative of the gravitational potential, 
t is the lookback time, with t = at the present and t — 
at the last scattering surface. The angular power spectrum 
of these temperature fluctuations (see Appendix A) is given 
by 
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where r is the comoving distance to lookback time, t, ji is 
the spherical Bessel function and P^^(k, r) is the 3-D power 
spectrum of $ fluctuations. To derive the final expression 
we have used Limber's approximation by assuming k f « l/r 
(|Limberll 19541 ; lKaiserlll992l ; IHu1|20001 ; IVerde et al.ll2000l . also 
see Appendix A) . 

The ISW effect consists of the temperature fluctuations 
described by these equations when linear theory is used to 
compute $ and its fluctuation power spectrum P§§- Using a 
simulation to determine the non-linear contributions we can 
quantify the full ISW plus Rees-Sciama effect. In Fourier 
space, the time derivative of the gravitational potential can 
be expressed as: 
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where a is the expansion factor, Ho is the Hubble constant, 
f2 m is the present mass density parameter and S is the time 
derivative of the density fluctuation. Combining this with 
the Fourier space form of the continuity equation, S(k,t) + 
ik ■ p(k, t) = gives: 
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where p(k,t) — [1 + S(k,t)]v(k,t) is the momentum density 
field in Fourier space divided by the mean mass density. This 
enables us to estimate the Fourier transform of the $ field 
of the simulation from the Fourier transforms of the den- 
sity and momentum fields. Using equation (|3} , the resulting 
power spectrum, P^(k,t) = (27r) _3 ($(fc, t)$*(k, t)), can be 
written as 
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Figure 1. Scaled $ power spectrum, P^^(k), at different red- 
shifts, A 2 (fc) = k 3 V^(k) /2tt 2 . The dotted lines are the mea- 
surements from the L-BASICC simulation. The solid lines are 
our model, while the dashed lines are the linear theory. The de- 
viation of the simulation results from linear theory happens at 
larger scales as redshift increases. We also find that the deviation 
from linear theory for the $ field occurs at larger scales than that 
of the density field at all rcdshifts. 



Figure 2. ISW angular power spectrum coming from different 
redshift intervals, where l(l + l)Ci /2n = (AT/T) 2 . The solid lines 
are given by Eq. J2j and Eq. J5}, evaluated using our model of the 
measurements from the L-BASICC simulation. The dashed lines 
are linear theory. The power spectrum at < z < 6 shows that 
the deviation of the simulation results from linear theory starts 
at I < 100. The deviation starts at smaller I as redshift increases. 



In linear theory, P ss (k,t) = k 2 P vv (k,t) = D(t) 2 Plj(k) 
and P ss (k,t) = kP Sv (k,t) = D{t)D(t)P ss n (k), where P 55 n (k) 
is the linear density power spectrum at the present time and 
D(t) is the growth factor normalised to be unity at present. 
Therefore, the power spectrum of the linear ISW effect is 



Ho 
k 



H(t)D(t)(l-(3) 



P)?s n (k), 



(6) 

where (5 = A d l ™ n ® m f^Si 6 (*) - For easy comparison at differ- 
ent redshifts in the simulation, we defined a scaled $ power 



spectrum, 
simply 



n m ] which from Eq. (5) is 
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Our measurements of the V^{k,z) power spectrum are 
shown in Fig. [T] The results from linear theory are also plot- 
ted. We find the total scaled "I? power spectrum can be well 
fitted by a broken power law plus the linear scaled 3> power 
spectrum = P|| Jin + P$, where 



nt in (k,z) = 



A 



(10fc)- 4 /0.75 + (lOfc)- 4 /^ 



(8) 



Here A and B are two free parameters that we use to fit 
the model to the simulation results at each redshift up to 
z = 6. To interpolate the model to intermediate redshifts 
we linearly interpolate the values of A and B from the near- 



est two simulation outputs. Our model is compared to the 
simulation results in Fig. [T] 

We see in Fig. [1] that the linear theory reproduces the 
ISW+RS P 44 at z = only at k < 0.1 h Mpc -1 . It fails 
at progressively larger scales as the redshift increases. By 
z = 2, linear theory agrees with the simulation results only 
at k < 0.02 h Mpc^ 1 . The reason for this surprising be- 
haviour is that the linear part of the drops quickly to 
zero as the relative importance of Qa diminishes at high 
redshift, while the non-linear part evolves more slowly with 
redshift. Therefore, the deviation of the total power spec- 
trum from linear theory happens at larger scales at higher 
redshifts. We find that the momentum power spectrum, P ss , 
and the correlation power spectrum of the density and mo- 
mentum, P sg , behave similarly to the P 44 power spectrum, 
namely, their deviation from linear theory occurs at larger 
scales at higher redshift. This is in contrast with the power 
spectrum of the density field which deviates from linear the- 
ory on progressively larger scales at lower and lower redshift. 
In another words, at the same redshift, the deviation from 
linear theory occurs at smaller scales for the density field 
than for the other fields. 

The sharp increase of measured from the simula- 
tion at small scales (k > 1 h Mpc -1 ) is due to discrete- 
ness in the 448 3 particle L-BASICC simulation. We used 
the m uch higher resolution 2160 3 particle Millennium simu- 
lation (|Springel et al.ll2005t ) to verify that our model remains 
accurate at smaller scales and is robust to shot noise correc- 
tions. 
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Figure 3. Evolution of the <I> power spectrum A 2 (fc) = k^P^^(k) /2tt 2 for specific spatial (left) and angular (right) modes. The solid lines 
show our model of the simulation results. The dashed lines show the results of linear theory. The power A 2 (fc) decreases monotonically 
as a function of z at all scales in linear theory. However, the total power seems to be independent of z at high redshift. The deviation 
from linear theory increases with z and k. The angular power spectrum A 2 (I) shows no deviation from linear theory up to I = 10. For 
I > 50, deviations appear at all redshifts. 



We can now compute the induced angular power spec- 
trum of CMB temperature fluctuations by performing the 
integral in equation (2) over the redshift range < z < 6 
using our model for the 3-D power spectrum, P^,^(k, z). The 
overall result is shown in Fig. [2] along with the contributions 
coming from different redshift intervals. For the overall an- 
gular power spectrum the deviation of the model from the 
linear theor y happens at I ~100. This result confirms the 
prediction of lCooravl i|2002bT ) based on the halo model. How- 
ever, we also see that the failure of linear theory, as judged 
by our simulation results, occurs at smaller and smaller / as 
redshift increases. For example, above z = 5, the deviation 
occurs at I < 20 and, for larger values of I than this, linear 
theory becomes extremely inaccurate. 



In order to evaluate how the breakdown of linear theory 
depends on redshift, we plot the evolution of the <3? power at 
a given scale as a function of redshift in Fig. [3] Generally, 
the deviations of linear theory from the simulation results 
decrease with scale and increase with redshift. At k — 0.01 h 
Mpc -1 , deviations start to be seen at z ~ 3 and, at k — 0.1 h 
Mpc -1 , linear theory has become inaccurate at all redshifts. 
In the right-hand panel, which shows results in / space, we 
find no deviations up to I ~ 10, but for I > 50, linear theory 
has clearly broken down at all redshifts. Interestingly, at 
high redshift, the $ power in the simulation appears to be 
independent of z while, in linear theory, this quantity drops 
monotonically with z. 



3 THE LSS-CMB CROSS-CORRELATION 

To illustrate the contribution of the Rees-Sciama effect 
to the cross-correlation of the density field S with <3>, 
we can compute the 3-D cross-correlation power spectrum 
P^ s (k, z) —< <fr(k, z)8* (k, z) > from our simulations: 

V^(k,z)=P ss {k,z)-^^- (9) 

In linear theory, V is (k, z) = D 2 (l~P)P l s T(k), where P l S s n (k) 
is the linear density power spectrum at the present time. Re- 
sults from our simulations are shown in Fig. [4] Comparing 
with linear theory, we find that the non-linear contribution 
appears at somewhat smaller scales than that of the auto- 
correlation power spectrum of At z — 0, the deviation oc- 
curs at k ~ 1 h Mpc -1 . However, it then dominates rapidly, 
making the cross-correlation power spectrum negative. This 
indicates that once the non-linear effect dominates, the po- 
tential of overdense regions evolves faster than the expansion 
of the universe, becoming deeper and thus imparting a net 
redshift to CMB photons passing through them. The sense 
of the effect from underdense regions is reversed, but the 
effect generated by the overdense regions is dominant and 
induces a negative cross-correlation between the CMB and 
the LSS.0 

3 An alternative way of seeing the sign of the Rees-Sciama effect 
is to consider the cross power spectrum P± s which from Eq. (9) 
can be seen to be proportional to d(Psg/a?)/dt. The ISW effect 
is the result of the linear density power spectrum growing less 
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Figure 4. Left: The cross power spectrum of <j> with <5 at different redshifts. The dotted lines represent the measurements from the 
L-BASICC and the dashed lines linear theory. The solid lines are our model fit to the sum of the linear theory and the non-linear 
contribution in the L-BASICC simulation. The non-linear effect begins to appear at k ~ 1 h Mpc" 1 at z = 0, and at larger scales at 
higher redshifts. It rapidly suppresses the cross power spectrum and causes it to become negative on small scales. Right: the same as 
the left panel but on a linear scale. 



To quantify the effect on the large-scale ISW cross- 
correlation measurements, we model the 3-D cross power 
spectrum from the simulations at each redshift output. We 
use linear theory to model the linear regime but once the 
cross power spectrum starts to deviate from linear theory, we 
fit it with a function of the form k 3 P^ s (k) = Ai+A2k+Ask 2 , 
where A\ , A2 and A3 are free parameters. Interpolating to 
intermediate redshifts, we are then able to calculate the pro- 
jected 2-D power spectrum. 

The cross-correlation between LSS and CMB maps has 
been shown to be a powerful tool for verifying the existence 
of dark energy and constraining its properties. Current mea- 
surements of the cross-correlation have low statistical sig- 
nificance because the volumes probed by LSS surveys are 
relatively small, but this situation will improve greatly with 
upcoming surveys. For example, Pan-STARRSl will survey 
three quarters of the sky, obtaining photometry for galaxies 
up to ~ 24.6 mag in the gr-band. The mean galaxy red- 
shift in this "3tt survey" will be z ~ 0.5. Pan-STARRSl 
will also carry out a deeper but smaller "MDS" survey cov- 
ering 84 sq deg of the sky to ~ 27.3 mag in g for which 



rapidly than a 2 once curvature or A become dominant. In con- 
trast the Rees-Sciama effect results from the non-linear growth 
of the density power spectrum, which for the times and scales o f 
interest grows faster than linear theory fe.g. ISmith et al1l2003t) . 
Consequently the contribution to the cross power spectrum from 
the Rees-Sciama effect has the opposite sign to that of the ISW 
effect. 



z ~ 0.8 l|Cai et al.ll2008l ). Cross-correlating such photomet- 
ric redshift galaxy samples with a CMB map (from WMAP 
or Planck) will make it possible for the first time to perform 
ISW tomography. Galaxy samples would be divided into dif- 
ferent redshift slices and each one cross-correlated with the 
CMB map. Values of the dark energy equation of state pa- 
rameter, w, could then be measured using the results from 
the different redshift slices, effectively constraining the evo- 
lution of w. 

To illustrate how ISW tomography may work, we follow 
iBaugh fc Efstathioul l|l995t ) and model the redshift distribu- 
tion of galaxies tracing the LSS as 

fO-z c ) 2 exp L(£^£cW2] if 2<2c 
N(z) (x l V ' P L v z o > J - ( 10 ) 

[0 if z > z c , 

but then choose the parameters zo and z c to emulate plau- 
sible photometric r edshift slices . (The same functional form 
was also taken by ICabre et all l|2006h to model the SDSS 
LRG sample.) We assume zo = 0.2 so that the width of N(z) 
is much greater than the expected photometric redshift er- 
rors. We shift the function into different redshift intervals 
by using z c = 0,1,2 and 3. The median redshift of these 
samples is z ~ 1.4zo + z c = 0.28, 1.28, 2.28 and 3.28 respec- 
tively. The cross-correlation power spectrum (derived in an 
analogous way to the auto-correlation function detailed in 
the appendix) is given as: 

C ?~° « 3 f L P ^ k = l -iz)Kz)N{z)H(z)/r 2 dz, (11) 
c Jo ' 

where P^ s (k, z) is the cross power spectrum of the potential 
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field and the galaxy density field, b(z) is the galaxy bias 
parameter at redshift z, and N(z) is the normalised galaxy 
selection function, where J N(z)dz = 1. We adopt the small 
angle approximation in which k — l/r(z), where r(z) is the 
comoving distance. For simplicity, in this illustration, we 
assume the galaxy bias parameter to be unity. In angular 
space, the cross-correlation becomes: 



,*-s 



(0) 



(12) 



where Pi are Legendre polynomials. In actual measurements 
of CMB fluctuations, the monopole and dipole are sub- 
tracted. Therefore, we set the power at I = and I — 1 
to zero before converting the signal into real space. To en- 
sure that the results at smaller angles (9 < 1°) converge 
accurately, we sum the power up to I = 10 000. 

The cross-correlation results are shown in Fig. The 
contribution from the non-linear RS effect can be seen to be- 
come increasingly important as the redshift of the sample in- 
creases. The cross-correlation power spectrum decreases and 
deviates from linear theory rapidly at I ~ 500 due to the non- 
linear effect. It turns negative a t I ~ 1000. This result i s con- 
sistent with that obtained by iNishizawa et all (|2008l ) who 
used a similar method to illustrate the impact of the ISW 
and RS effects on the CMB-weak lensing cross-correlation. 
In angular coordinates, shown on the right-hand panel, the 
RS effect is negligible at 9 > 1° at all redshifts. For the high 
redshift samples, it becomes important at sub-degree scales 
for the high redshift samples where it suppresses the cross- 
correlation power spectrum by about 5 — 10% at acrminute 
scales. 

The statistical significance of current measurements of 
the CMB-LSS cross-correlation is not yet high enough to 
detect the effects we are discussing. Most current measure- 
ments can only det ermine the cross-c orrelation at degree 
scales or above ("e.g. ICabre et al.ll2006l ). Future CMB or SZ 
surveys with high resolution might be able to resolve this 
contribution. 



the importance of the Rees-Sciama effect may become sig- 
nificant for high redshift samples. The effect of the non- 
linear cross-correlation at scales of arcminutes would con- 
taminate the SZ signal in the CMB and thi s could con- 
fuse i ts interpretation (e .g. Fosalba et al.l|2003l; Diego et "all 
| 2003l; iMvers et alj 12004 : iLieu et al.l l20od : ICao et all 120061 ; 
bielbv fc Shanksll2007h 

Our analysis is based on a simulation that assumes a 
ACDM cosmology. The non-linear contribution depends on 
the values of the cosmological parameters. In a flat uni- 
verse with a cosmological constant, the RS effect will be- 
come increasingly dominant relative to the ISW effect as 
the value of Qa decreases. In the most extreme case, if 
SIa = 0.0, the ISW effect will vanish, leaving only the RS 
effect | |Seliaklll996h . The analysis of this paper could be gen- 
eralized either using re-norm alised perturbation theory (e.g. 
ICrocce fc Scoccimarrol [20061 ) . or simulations with different 
dark energy models. In any case, more general modelling 
of the non-linear effect will be required for an accurate in- 
terpretation of future measurements of the LSS-CMB cross- 
correlation. 
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4 CONCLUSIONS 

We have used an N-body simulation to calculate the non- 
linear (Rees-Sciama) contribution to the Integrated Sachs 
Wolfe effect. The comparison of the 3-D and 2-D power 
spectra measured from the simulation with those given by 
linear theory reveals a strong nonlinear contribution whose 
physical scale increase with redshift. We investigated the 
strength of this effect on the cross-correlation of the CMB 
with galaxy samples in terms of angular power spectra and 
in angular coordinates at different redshifts. We find that 
there is a non-linear contribution to the cross-correlation sig- 
nal at sub-degree scales. The non- linear effect alters not only 
the amplitude, but also the shape of the cross-correlation 
power spectrum. With current galaxy samples which cover 
relatively small volumes, it is not yet possible to disentan- 
gle the contribution of the RS effect from that of the ISW 
effect within the noise. However, in future surveys like Pan- 
STARRS and LSST, for which the number of galaxies and 
the sky coverage will increase dramatically, the error bars 
on the cross-correlation will be much smaller. In this case, 
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Figure 5. Left: The cross-correlation power spectrum of galaxy samples at different z with the CMB. The dashed lines are given by 
linear theory. The solid lines are the sum of the linear theory and the non-linear contribution, which is given by our model fitted to the 
L-BASICC simulation. The non-linear effect begins to appear at / ~ 500. It rapidly makes the cross power spectrum become negative. 
Right: The cross-correlation of galaxy samples at different z with the CMB in angular space. The non-linear effect suppresses the ISW 
effect at sub-degree scales. It is negligible at low redshifts and large scales. 
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APPENDIX A: ANGULAR POWER SPECTRA 

Here we derive the relationship between the 3-D power spec- 
trum of gravitational potential fluctuations, P^^(k,t) = 
(27r)~ 3 {"l>(fc, £)$*(£;,£)) and the resulting angular power 
spectrum of the induced CMB temperature fluctuations. Ex- 
panding the pattern of temperature fluctuations, AT(r)/To, 



8 Cai et al 



in terms of spherical harmonics we have 
AT(f) 



aim = 



T 



Yim(r)dr 



which using equation |T} becomes 



(llr 



~# I Y ^f) / ${r,t)dtdr. 



(Al) 



(A2) 



Writing <£>(f , t) in terms of a Fourier expansion and using the 
spherical harmonic expansion of a plane wave, exp(ik ■ r) = 
47r I]i m *' 1 3i{kr)Y lm (k)Yi m (f) this becomes 



dir. 



(2tt) 3 c 2 
2 x 4tt . 

ns3 / y '-W 



Ylmif) I I <$>{k,t)exp(ik ■ r)dkdtdr 



fry* J ' mv 'j tt 7 * (M)x 

2J iji {kr)Yi* m , (fc)Yi' m ' (f) dfc dt df 



(see al so lLimberlll954l; iKaiserl Il992l ; IHuI |2000| ; IVerde et all 
1200(f) . IVerde et al l <|2000h find that the difference between 
this approximation and the full calculation is less than 3% 
at / > 20. In this paper we are mainly concerned with even 
smaller scales (I > 500) and so we are justified in using 
Limber's approximation. We also tested the difference be- 
tween using k = l/r and t he more accurate k = (I + l/2)/r 
l|Loverde fc Afshordillioog) and find very little difference to 
our results. 



t)i ji (kr)Y* m (k) dk dt. 



Hence the angular power spectrum, Ci, is given by 



11' $mrr> 


' = {dirndl' 
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$(fc, r)i ji{kr)Y* m {k) dk dr 



^*{k' ,r')i 1 ' ii{k'r')Y llm ,{k')dk'dr' 
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Using the identity (&{k)$>* (k')J = (27r) 3 <5(fc - k')P^(k) 
and the orthogonality relationship of spherical harmonics 

/ Y*m{f)Yl'm'{f) df = 8 w 8 mm , (A5) 

J 4tt 

this becomes 

8 2 f f TL r rL 

ClS U lS mm l = J- / / P®i,{k,r,r')i jl{kr)Y,* m {k) 

x i 1 ' ii(kr')Y llm ,(k) dr' dr dk 
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k 2 P^{k,r, r')ji(kr)ji(kr')dr'drdk5u>6r, 



Q = 



— f r 

^ J Jo 



k 2 P^^(k, r, r')ji(kr)ji(kr')dr'drdk. 

(A6) 

This exact relationship can be simplified by using Limber's 
approximation. For small angular separations, 8, at comov- 
ing distance, r, the wave number, k, can be expressed in 
terms of its components parallel and perpendicular to the 



line of sight and approximated by k = + k'j_ ~ k±, 

where k± — 2n/r6 ps l/r 3> fey ~ 1/Ar, namely, the power 
is dominated by that perpendicular to the line of sight and 
there is no correlation between different shells of Ar along 
the line of sight. Combining this with the orthogonality re- 
lation for spherical Bessel functions, 



— / k 2 ji(kr)ji(kr')dk = Sir — r')/r' 2 
7T , 



we arrive at 
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